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ABSTRACT 

The temporal behaviour of segmentation clock oscillations show phase synchrony via mean 
field like coupling of delta protein restricting to nearest neighbours only, in a configuration 
of cells arranged in a regular three dimensional array. We found the increase of amplitudes 
of oscillating dynamical variables of the cells as the activation rate of delta-notch signaling is 
increased, however, the frequencies of oscillations are decreased correspondingly. Our results show 
the phase transition from desynchronized to synchronized behaviour by identifying three regimes, 
namely, desynchronized, transition and synchronized regimes supported by various qualitative and 
quantitative measurements. 

.ftTEYWORDS : Delta-notch coupling, phase synchrony, Hilbert transform. Phase locking 
value, mean field coupling. 

PACS numbers: 



I. INTRODUCTION 



The origin of oscillatory behaviour of segmentation 
clock genes is believed to be due to negative feedback 
genetic regulation of their own products [iHll- For ex- 
ample, in vertebrate developcment, periodic formation 
of somites during presomitic mesoderm is considered to 
be responsible to exhibit oscillatory expression by these 
genes during sometogenesis and their period of oscillation 
is very close to the period of segmentationQ-Q- It has 
been reported that during this signal processing neigh- 
bouring cells are in contact and these cells communicate 
among the neighbouring cells showing synchronization in 
the dynamical variables 0, • 

There have been various theoretical studies regard- 
ing the dynamical stability, periodicity and synchroniza- 
tion of the various oscillating variables of segmentation 
clock Since synchronization and oscillation in 

the clock gene regulations are believed to be responsi- 
ble for somite boundary, these properties are considered 
to be essential in segmentation process [l^. If a group 
of such cells are considered, signal processing among the 
neighbouring cells are done via delta-notch signaling pro- 
cess to achieve synchronization among them |16l . llSl ll9[ . 
However, since the delta proteins cannot diffuse from one 
cell to another, the incorporation of delta-notch coupling 
in the theoretical models so far developed is still ques- 
tionable. One attempt in this direction without consid- 



ering time delay was to introduce a coupling term in the 
model associated with activation rate of delta-notch sig- 
naling which enable to control the fluctuation in delta 
protein concentration in its neibouring cells [l^, [2lj . It 
has been reported that the cells in this coupling mecha- 
nism exhibit phase synchrony in one and two dimensional 
array of such cells |19l - [2l| . 

In this work we focuss on some issues namely, measure 
of phase synchrony among the cells and generalization of 
coupling mechanism into some more real situation. In 
the section II, we study the model of segmentation clock 
gene expression due to Uriu et.al. (2009). We generalized 
the model in three dimensional situation with coupling 
mechanism. Then we explain various methods to identify 
phase synchrony among the coupled cells. The results 
and discussions of our simulation results are presented 
in section HI. We draw some conclusions based on our 
results in section IV. 



II. MATERIALS AND METHODS 

In vertebrate sometogenesis, the expression of her gene 
in single cell model with and without time delay is es- 
sential in order to generate oscillatory behavior of the 
segmentation clock genes [l^, [2^ [2^ [2^ . Here, we con- 
sider a model in which simple kinetics of mRNA, protein 
in cytoplasm and protein in nucleus are considered with- 
out taking into account any time delay [2l[. The detail 
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FIG. 1: The schematic diagram of reaction network of cell to 
cell interaction via delta-notch coupling mechanism. 
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FIG. 2: The three dimensional array of cells showing nearest 
neighbour interaction. 



molecular reaction pathways are shown in Fig. 1. 

In this model [2l|, if Af, Y , Z and W indicate con- 
centrations of Her mRNA, Her protein in cytoplasm, 
Her protein in nucleus and Delta protein respectively in 
a single cell, then the time evolution of these variables 
without considering time delay in describing the kinetics 
in the cell are given by. 
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the differential equations arc given as follows: vi is the 
basal transcription rate, Vc indicates the activation rate 
due to Delta-Notch signaling, V2 is maximum degrada- 
tion rate of Her mRNA, refers to translation rate of 




FIG. 3: The three dimensional plot of M{t), Y{t) and Z{t) for 
Vc = 0.263, 0.5 and 0.7. The other values of the parameters 
used in the simulation are; ^1 = 0.0135, V2 = 0.6, 1^3 = 0.575, 
VA = 0.851, i/5 = 0.021, = 0.162, = 2.465, vs = 9.583, 
fci = 0.157, fca = 0.104, = 0.142, fcg = 0.13, kr = 0.49, 
fcg = 9.72,n = /i = 2 [3. 



Her protein, 1/4 is the maximum degradation rate of Her 
protein in cytoplasm, 1/5 indicates transportation rate of 
Her protein, is maximum degradation rate of Her pro- 
tein in nucleus, is the sysnthcsis rate of Delta protein, 

is the maximum degradation rate of Delta protein, ki 
is threshold constant of supression of Her mRNA tran- 
scripted by Her protein, k2 is Michaelis constant for Her 
mRNA degradation, ^4 is Michaelis constant for Her pro- 
tein degradation in cytoplasm, fcg is Michaelis constant 
for Her protein degradation in nucleus, k-j is threshold 
constant for the suppression of Delta protein synthesis by 
Her protein, fcg is Michaelis constant for Delta protein 
degradation in nucleus. 

A group of such cells do communicate via delta-notch 
coupling mechanism when the cells are very close to 
each other or almost in contact each other. Since Delta- 
proteins are confined only on the cell membrane and can- 
not diffuse from one cell to another, signal is considered 
to pass from one cell to another by activating the delta- 
notch signaling pathway when the cells are in contact or 
very close. We consider [N x N x N) cells which are ar- 
ranged in a regular three dimensional array in this type 
of communication, assuming the sizes of the cells are the 
same as shown in Fig. 2. In this situation, any cell is 
allowed to couple with six nearest neighbour cells. Since 
the coupling variable W is common to (N x N x N) 
identical cells, this species can be an effective means of 
coupling each cell to every other via i X^i^i where, 
L — 6 and summation is restricted to nearest neighbours 
only for every cells. The common species provides a "lo- 
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FIG. 4: Plots of the variables M{t), Y{t), Z{t) as a function 
of time. The same parameters are used as in Fig. 3. 



calizcd" mean field like coupling whereby the cells com- 
municate, and the remaining variables synchronize. 

For one and two dimensional array, the allowed num- 
ber of coupled cells are two and four respectively. So 
for the cells in this configuration, the delta-notch signal- 
ing activation rate, Vc plays switching on or off of signal 
processing among the cells and in "on switch" condition, 
coupling in any cell is done via coupling of its nearest 
neighbour cells. Therefore in delta-notch signaling pro- 
cess, we have the following coupled oscillators. 
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where, a,b,c ^ 1,2, . . . , N , L ~ 2, 4, 6 for 1, 2, 3 dimen- 
sional array of cells and (a, a), {a, /3), (c, 7} indicate that 
the pairs (a, a), (6, ,9) and (0,7) respectively are near- 
est neighbours. The mode of coupling among the cells 
is quite similar to mean field coupling scheme (2^ . [25j 
but the way how the coupling is incorporated among the 
cells in this case is localized within the nearest neighbours 
only. In this coupling mechanism, once the coupling is 
switch on with a particular value of Vc, then the cells 
start synchronized among themselves to exhibit coordi- 
nated behaviour. In our simulation we use fixed bound- 



ary conditions in the coupling, i.e. Wobe ~ 0, WaOe = 0, 

WabO = 0; Wn+IM = 0, Wa,N+l,e = 0, Wab.N+1 = 0. 

It has been pointed out that the identification of phase 
synchrony of any two identical systems can be done qual- 
itatively by the measurement of the time evolution of in- 
stantaneous phase difference of the two systems [2^-[2^. 
It is possible if one defines an instantaneous phase of an 
arbitrary signal x{t) via Hilbert transform j27| 
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where P.V. denotes the Cauchy principal value. Then, 
one can determine an instantaneous "phase" (j){t) and 
an instantaneous "amplitude" A{t) of the given signal 
through the relation, A(t)e''^(*) = x{t)+ix{t). Given two 
signals of two systems Xi (t) and Xj {t) , one can therefore 
obtain instantaneous phases (t>i{t) and (j^jit)] phase syn- 
chronization is then the condition that A(f>ij = m(f>i — n(f>j 
is constant, where m and n are integers. Of most interest 
are the cases A(j)ij = or tt, namely the cases of in-phase 
or anti-phase, but other temporal arrangements may also 
occur. 

The phase synchronization of the two identical systems 
can also be identified by doing pecora-caroll type recur- 
rence plot of the variable of one system, say x with the 
corresponding variable x' of the other system simulta- 
neously on two dimensional cartesian co-ordinate system 
[30| . The rate of synchronization between the two sys- 
tems can be determined by the rate of concentration of 
the points in the plot along the diagonal. If the two 
systems are uncoupled, then the points on the plot will 
scatter away randomly from the diagonal. 

The rate of phase synchrony of the two systems can be 
estimated quantitatively by measuring the "phase lock- 
ing value" (PLV) of the two signals of the two systems 
3lj . The phase locking value, which is used to quan- 
tify the degree of synchrony, characterizes the stability 
of phase differences between the phases (piit) and (/>j(t) 
of two signals Xi{t) and Xj{t) of itli and jth systems and 
can be defined within a time period T by. 



Pit) = ^ 



t-T 



(10) 



The range of PLV value is [0,1]. When the value of PLV 
is zero, the two systems are uncoupled, whereas if the 
value of PLV is one then the two systems are perfectly 
phase synchronized [sij . 



III. RESULTS AND DISCUSSIONS 

We first present one time simulation results of M, y 
and Z of a single cell as a function of time for the values 
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FIG. 5: Plots of time evolution of M{t), Y{t), Z{t) of 10 oscillators out of (10 x 10 x 10) three dimensional array of oscillators 
are shown in (a), (c) and (e). The corresponding plots of A(j}ij{t) vs time are shown in (b), (d) and (f) respectively. 



of i^c = 0.263, 0.5 and 0.7 respectively in Fig. 3 by solving 
differential equations 1,2,3 and 4 by standard 4th order 
Runge Kutta numerical integration method (32l |. In the 
3D plot, the dynamics shows three different limit cycle 
oscillations for these three different i/c values. In Fig. 4, 
the corresponding dynamics of the three variables for two 
coupling cells are shown in panels (a), (c) and (e) indicat- 
ing phase synchrony irrespective of variation in i^c- The 
claim is supported by the phase plots i.e. A0i2 vs time 
in Fig. 4 (b), (d) and (f), along with pecora-caroll type 
recurrence plots in the panels (g), (h) and (i). It is no- 
ticed from the plot that the amplitudes of oscillations of 
all variables are increased as Vc is increased however the 
frequencies of oscillations are decreased correspondingly. 
The phase locking value for whole time period (0-1500) 
minutes using equation (fTO|) is found to be 0.97±0.015 



indicating fairly well synchronized. There are signatures 
of fluctuations in the beginning of each plot due to the 
cells evolve in different initial conditions and takes some 
time to synchronize. This means that once the coupling 
is switch on via Vc among the cells, the variation in Vc 
will hardly effect in the synchrony rate of the cells but 
will occur at different amplitudes. 

Next we take (10 x 10 x 10) such cells arranged in three 
dimensional array. We take i^c ~ 0.263 and apply fixed 
boundary conditions and simulate differential equations 
5, 6, 7, 8 and the results of M, Y and Z function 
of time for 10 cells are shown in Fig. 5. The coupling 
variable for each cell is W via mean field like coupling 
mechanism restricting to nearest neighbours only and the 
coupling is switch on at time 300 minutes. In the plots of 
Fig. 5 (a), (c) and (e), three regimes i.e. desynchronized. 
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FIG. 6: The pecora-caroU type recurrence plots of the M, Y and Z of 10 oscillators shown in Fig. 3 are shown in the upper 
three panels. i,j = 1, 2, 10, i ^ j. The corresponding P values of the three variables for different time periods are shown in 
the lower three panels. 



transition and synchronized regimes are seen clearly. The 
phase transition like behaviour is supported by the phase 
diagrams (A0ij vs time) of pairs of cells in figures (b), 
(d) and (f) where desynchronized regime is indicated by 
random evolution of A.(j)ij with time, transition regime is 
identified by weak fluctuation of around a constant 
value and synchronized regime is indicated by constant 
A(f>ij with time. Again this claim of phase transition 
is verified qualitatively by pecora-caroU recurrence plots 
in the upper three panels of Fig. 6. In these plots if 
the cells are desynchronized, the points in the plots will 
spread away from the diagonal, however in the transition 
regime the points start concentrating towards the diago- 
nal. When the points lie along the diagonal the cells are 
synchronized as shown in the plots. 

Fixing the configuration of cells, we calculated A(/)ij 
of pair of far distant cells (i.e. i=l, j=10) func- 
tion of time, which are already included in the respective 
diagrams, and found the same pattern of phase transi- 
tion. This means that once the coupling is switch on. 



all the cells process information almost instantaneously 
even at far distant cells also in the configuration. Hence 
long range information transfer or relay is done instanta- 
neously when the cells are coupled. 

The quantitative identification of this phase transition 
can be explained more clearly by the measurement of 
phase locking values {Pm, Py and Pz) of the variables 
A/, Y and Z for different time intervals by using the equa- 
tion ([T0|) . The results are shown in the lower three panels 
of Fig. 6. In these plots, it is seen distinctly the demar- 
cation of the three regimes, desynchronized, transition 
and synchronized regimes. In the synchronized regime, 
the values of Pa/, Py and Pz are 0.994±0.002 starting 
from time interval (400-500) minutes to (1100-1200) min- 
utes. In the time intervals before switching the coupling, 
the values of Pa/, Py and Pz are in 0.77±0.05 indicating 
desynchronized regime. But in the time intervals 300 to 
500 minutes, the values of Pa/, Py and Pz monotonically 
increased indicating transition regime. 
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IV. CONCLUSION 

The phase synchrony shown by temporal oscillations of 
the variables indicates that once the coupling is switched 
on with the mean field like coupling mechanism with a 
value of activation rate above some critical value to 
exhibit synchrony, then the increase of Vc does not affect 
much in the rate of phase synchrony. However it affects 
in the magnitudes of amplitudes and frequencies of oscil- 
lations of the dynamical variables of the cells. It is also 
to be noted that if the variables evolve with time from 
different initial conditions, after coupling is switched on, 
it takes some time to get synchronized. 

Since the cells coupled only when they arc very close 
or in contact and the delta protein cannot diffuse from 
one cell to another, the mean field like coupling scheme 
is the most reasonable way of signal processing in this 
situation. We take three dimensional array of cells to 



see the behaviour of phase synchrony in more real situa- 
tion. Our results interestingly show phase transition like 
behaviour by identifying desynchronized, transition and 
synchronized regimes. Relay information transfer among 
the far cells is done among the cells themselves once the 
coupling is introduced in a group of cells. To investigate 
it in more real situation one has to incorporate noise 
fluctuations arising from internal molecular events and 
external environmental fluctuations. 
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